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Abstract 


An unavoidable step in the process of space exploration is to use high-power, very large 
spacecraft launched into Earth orbit. Obviously, the spacecraft will need powerful energy 
sources. Previous experience has shown that electrical discharges occur on the surfaces of 
a high-voltage array, and these discharges (arcs) are undesirable in many respects. 
Moreover, any high voltage conductor will interact with the surrounding plasma, and that 
interaction may result in electrical discharges between the conductor and plasma (or 
between two conductors with different potentials, for example, during docking and 
extravehicular activity). One very important aspect is the generation of electromagnetic 
radiation by arcing. To prevent the negative influence of electromagnetic noise on the 
operation of spacecraft systems, it seems necessary to determine the spectra and absolute 
levels of the radiation, and to determine limitations on the solar array bias voltage that 
depend on the parameters of LEO plasma and the technical requirements of the spacecraft 
equipment. This report describes the results of an experimental study and computer 
simulation of the electromagnetic radiation generated by arcing on spacecraft surfaces. A 
large set of high quality data was obtained during the Solar Array Module Plasma 
Interaction Experiment (SAMPEE, flight STS-62) and ground test. These data include the 
amplitudes of current, pulse forms, duration of each arc, and spectra of plasma waves. 

A theoretical explanation of the observed features is presented in this report too. The 
elaborated model allows us to determine the parameters of the electromagnetic noise for 
different frequency ranges, distances from the arcing site, and distinct kinds of plasma 
waves. 


’National Research Council-NASA Research Associate at Lewis Research Center. 
*Chief, Space Environment Effects Branch. 

^Physicist, Space Environment Effects Branch. 
s Professor, Tuskegee University, Tuskegee, Alabama. 
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Nomenclature 
a arc length, m 

B magnetic field strength, Gs 

C capacitance, F 

E, E electrical field strength, V/m 

/ frequency of plasma density fluctuations, Hz 

fm lower hybrid resonance frequency, Hz 

fm upper hybrid resonance frequency, Hz 

F e electron gyrofrequency 

F t ion gyrofrequency, Hz 

I, i arc current, A 

Imax arc peak current, A 

J, J electric current density, A/m 2 

k, k wave vector, m* 1 

ks Boltzman constant, J/K 

n e electron number density, m' 3 

w, ion number density, m' 3 

N arc number of arcs per one experiment 

p neutral gas pressure, Pa 

R distance from spacecraft surface, m 

1 time , s 

T e electron temperature, eV 

T g neutral gas temperature, K 

T t ion temperature, K 

U LP output voltage, V 

v e electron thermal speed, m/s 

V bias voltage, V 

y the decrement of plasma waves, s' 1 
At duration of an arc, s 

q> the angle between the wave vector k 

and the Earth magnetic field B 
0(0 arc current spectrum, A/MHz 

X the wave length, m 

Ac the Debye length, m 

«9 angle between the Earth magnetic field and the arc current 
a cross section, m 2 

texp duration of individual experiment with arcing, s 
o)=2jrf frequency, Hz 

1 . Introduction 


Forty years experience of space exploration show that large spacecraft, such as Space 
Station Freedom, will be deployed in low Earth orbit to provide a basis for further 
achievements in space science and space technology. Obviously, large and complicated 
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structures demand powerful energy sources to keep all systems in operating condition. 
We believe that high-voltage solar arrays will serve as the primary electrical energy 
sources for a long time in the future. However, previous experience has shown that 
electrical discharges occur on the surfaces of a high-voltage array M , and these discharges 
(arcs) are undesirable in many respects. One very important aspect is the generation of 
electromagnetic radiation (EMR) by arcing. 1 * * * 5 To prevent the negative influence of 
electromagnetic noise on the operation of spacecraft systems, it seems necessary to 
determine the spectra and absolute levels of the radiation, and to determine limitations on 
the solar array bias voltage that depend on the parameters of LEO plasma and the 
technical requirements of the spacecraft equipment. 6 Moreover, any high voltage 
conductor will interact with the surrounding plasma, and that interaction may result in 
electrical discharges between the conductor and plasma (or between two conductors with 
different potentials, for example, during docking and extravehicular activity). 

A number of experiments have been done during last decade to determine arcing rate and 
threshold, the magnitude and duration of current pulses, the dependence of arcing rate on 
bias voltage and plasma density, and the spectrum and intensity of generated 
electromagnetic radiation. An extensive review of space experiments is published in the 
paper by Hastings 7 , and new data obtained in the last two years may be found in papers 
by Hillard and Ferguson 8 and Kuninaka 9 . The advantage of space experiments is 
obvious: solar cells (or samples of solar cells) operate in the natural plasma environment 
during many days or even months, and both short time and long time processes caused by 
interaction with plasma can be analyzed . But there are always strong limitations on the 
volume and types of data that might be collected because of technical requirements for 
experimental devices employed on spacecraft. On the contrary, ground experiments have 
the advantage of a controlled environment and unlimited choice of measurements that may 
be done. In other words, these two kinds of experiments are complementary to each 
other, and we use the results of the Solar Array Module Plasma Interaction Experiment 
(SAMPIE) 10 and the Langmuir Probe Experiment (LPE) 11 to create and verify a 
theoretical model that allows us to determine the parameters of the electromagnetic noise 
for different frequency ranges, distance from the arcing site, and distinct kinds of plasma 
waves. 

2. Basic equations 

The discharges that occur between negatively biased solar cells and the surrounding 
plasma can generate electromagnetic radiation with a broad spectrum and considerable 
amplitude. In order to obtain an appropriate expression connecting the discharge current 
characteristics with EMR spectra we start from Maxwell’s equations: 

1 dB An 1 d& . 

curlE = curlB = — J + — ; divE= 4 -ttp; divB = 0 (1) 

O ut, O C ut 

If we consider waves of length much shorter than the scale of the plasma density gradient, 

and if we disregard changes in the electron number density due to arcing, we can 

apply the Fourier transformation to the equation set (1): 
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E(x,f) = — — 7T • f E(a>,k) • exp(-/<y/ + /'kx) • d^kda) 


{l7t) 


( 2 ) 


It is easy to obtain the following relation between the Fourier transformations of the 
electrical field strength and the current density: 


f 2 

£ a 

2 


E(<y,k) +k -(k • E(/y,k)) = _ j(a,k) 


(3) 


The equation (3) allows us to find (in principle) the electromagnetic field if the source 
current is known. However, the main problem is to determine the current density 
J(o>,k) for all of space. Obviously, we can represent the current density as the sum of two 
terms: i) arc current; ii) current that is induced in the plasma by the electric field. 

The spatial distribution of the current density has the following form 12 : 


J(x,t)=J,(x,t)H(x) + J P (x,t)[l- H(x)] 


(4) 


where H(x) is a step function defined by the expressions: 

H(x)=l for all points x inside arc, and H(x)=0 elsewhere. 

Within the framework of the linear theory, the current in the plasma is proportional to the 
field strength 

Jp(©,k) = o(<y,k)E(fij,k) (5) 

where <r(e>,k) = cr^.(<i),k); (/, j - 1,2,3) is the conductivity tensor . 

The Fourier transformation of the current density can be written in the form: 


j(e>, k) = — — =-• J //(x)[J a (x, /) - J D (x, t )] ■ exp (iat - /kx) d^xdt+ 
{In r 

- — ~2 • ! Jp (x / /) • exp (icot - /kx) • d^xdt 


( 6 ) 


The last term in Eq.(6) is equal to J p (co,k). The first term in the Eq. (6) can be simplified 
by taking into account the inequality J P «J» that is valid inside the arc. 

The dielectric tensor of the plasma and the conductivity tensor are related according to 
the equation : 13 

K(<»,k) = i+— -CT(cy,k) (7) 

0) 


5 



where I = Sy is the identity matrix. 
Finally, the master equation has the form: 



The dielectric tensor has a complicated structure even for a cold magnetized plasma. 14 It 
can be written in the following form: 


K ij = Pb i 


bj+S 


■fa-vjh 


iDe ijk b k 


( 9 ) 


where P, S, and D are functions of frequency only. These functions can be represented as 
expressions with a strong dependence on the plasma resonance frequencies: 12 


P=l- 



_ (^ 2 

[f 2 -F?)\f 2 - F e] ’ 


D = - 


ffe'Fe 


(/ 2 -^ 2 )*(/ 2 




s ijk * s com P letel y antisymmetric tensor of rank three, and /= (o/2n. 


It should be noted that tensor K in the equation (8) doesn’t include an Ion Acoustic Wave 
branch because the ion temperature T ioa =0 by definition of the cold plasma. 

The formal solution of the Eq.(8) can be obtained in general form: 



( 10 ) 


where T = — -k-k 2 -i + k •(*•), and T T _1 = I (11) 

c 

The dispersion relations for electromagnetic waves are determined by zeros of the 
determinant of the tensor Tjj , and the electromagnetic field can be calculated by using the 
inverse Fourier transformation of the expression (10). However, to perform these 
calculations much algebra must be done. It seems more reasonable to calculate the 
electromagnetic field by using the approximation of the dielectric tensor within different 
frequency ranges . But, before starting this procedure we derive the Fourier components 
of the arc current density. 

3. Arc current 

The determination of the right term in the Eq. (8) requires a physical model of the 
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arc current. As it was observed in some experiments by measurement of UV radiation 
from arcs 3 , the diameter of the discharge region is very small («500 pm), and this value is 
much less than any wavelength that we suggest to be considered (A>10 cm). Moreover, 
we can estimate the upper limit of the arc cross section as 


s„=— 


nev. 


The amplitude of current / ra depends on many factors but it has never been higher than 
~10A for numerous space and ground experiments. 10 ’ 11 If we suggest that neutral gas 
should be fully ionized in the arc channel, the electron number density for LEO flight 
conditions is calculated according to the equality : n, = N& 10" -10 1J cm" 3 . The 

actual number density of electrical charge carriers should be much higher because 
secondary ionization caused by collision of electrons with neutral atoms (molecules) seems 
possible if one takes into account the desorption of gas from the cell surfaces. 15,16 The 


r 2eV' m 


» cm 

* 10 * — 

s 


electron velocity is certainly less than the maximum value = 

v m . J 

Thus, we get the following estimate for an upper limit on the arc diameter: 
d « iS^ ^ 1 cm . Obviously, the real diameter is much less than the estimate above. 
Therefore, the arc current can be written in the following form: 


= <&(*) ■ J{z) • £(x) • £(>>)• n. 


( 12 ) 


where 8(x) is the Dirac delta function, and n z is the unit vector in the z-direction. 

The temporal evolution of the arc current was studied in numerous experiments. 3,11 ’ 17,18 
Discharges on negatively biased solar cell samples were investigated during the space 
experiment SAMPIE . 8,1 ° The results of the statistical analysis and real curves &(t) will be 
shown below, but to an order of magnitude we know that the peak current is ~ (1-10)A, 
and the pulse duration is ~ (1-20) ps. 

We model the function <b(f) as in Fig. 1 . 

The current pulse form can be characterized by three parameters: rise time r, plateau 
duration At, and exponential relaxation time T. Of course, the spectrum 0(f) depends 
on these parameters. Computer simulation shows that the spectrum falls by 40 dB at 
frequency f=A0 MHz if r=0, and the spectrum goes down by 60 dB at the same frequency 
if t=0. 1 ps. An example of the spectrum is shown in Fig.2. 

The wave form of a current pulse was measured with sampling interval 5 ns in the 
course of LP experiments. 11 One example of such measurements is shown inFig.3 
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Fig. 3. Current pulse wave form measured during the ground experiment LPA_7. 

It is seen that the rise time is approximately equal to 3 ps , and the pulse duration slightly 
exceeds 10 ps. These data can be explained by the special experimental setup when an 
additional capacitor was installed between the anodized aluminium plate and ground to 
simulate the capacitance of real large solar array. The Fourier spectrum of current pulse is 
demonstrated in Fig. 4. The essential difference between the model spectrum (Fig. 2) and 
observed spectrum (Fig. 4) lies in the high-frequency region. The magnitude of harmonics 
does not fall below 70 db for frequencies f< SO MHz. It is caused by short time 
fluctuations in the arc current, and it is important for further computations because the 
plasma frequency does not exceed 20 MHz for LEO conditions. 

If we suggest that the arc channel has a simple structure like straight line perpendicular to 
the conductive plate , the spatial Fourier transformation of the current density can be 
represented by a simple formula: 19 

J(k) = j-(sin£ 3 a + /-(l-cosA: J a)) (13) 

*3 

where a is the arc length. 

It is easy to see that the current J(k) has maxima at wave numbers kiw=(2m+ l)/a 
(m-0,\,2...). We may conclude that an arc operates like an ordinary half-wave vibrator 
immersed in the plasma when we consider the radiation of elecromagnetic waves. 
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Fig. 4. The spectrum of current pulse shown in Fig. 3 


4. High frequency transverse waves 


We now begin our analysis of the electromagnetic radiation generated by arcing in the 
form of high frequency transverse waves. For LEO plasma, this means waves with 

frequencies f»f e &20 MHz and wavelength limited by the value X= x «15/». 


Because the length of an arc is not more than 10 cm*, the consideration below is valid 
within the frequency range 20</<3000 MHz. The frequency spectrum of the electric 
field strength can be obtained from equation (10) after performing the inverse Fourier 
transformation: 


E(x,<y) 



0(g)) - a 



• exp(/kx) 


(14) 


*This estimate is based on the experimental data obtained by Leung 17 and our calculations of the arc 
length. According to Ref.(6) , the arc length might be about equal to the plasma sheath size. This 
means that the parameter a could be as large as one meter if the bias voltage reaches several hundred 
volts. We don’t have the final answer in this problem. 
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If we define the angle 0 = cos 




form 


we may write the strength of electric field in the 


4n®(a)) -co - a 


£(W, ©,<*>)= w<c 2 


1 - cos' 


0+(3cos 2 0-l)-^ 


ic 


to 1 • hcl 2 a- [x|J 


(15) 


To compare the calculated value (15) with available experimental data and the technical 
requirements on the level of the electromagnetic noise, it is convenient to represent the 
field strength in the units of dB*pV/m*MHz. It can be written in the following form: 

E(x,<y) = 24O + 2O-log f £ ^ >0, ^l (16) 

v 4 ne n J 


where eo=8.85*10' 12 Cl/V*m is the dielectric constant. 

The spectrum (16) is shown in Fig. 5. To draw this chart we suggest the angle 0=x/2, the 
current amplitude I nux =l A, the arc length a= 10 cm, and the distance x=l m. As a result, 
we have a quite good agreement with experimental data obtained by Leung s . 
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Fig. 5 The frequency spectrum of the electric field strength at distance 1 m from 
an arc with the peak current 1 A. □ - experimental data from Ref.(5). 


It should be noted that all of the oscillations that are seen in Fig. 5 were not measured in 
the mentioned experiments because the highest values of amplitudes from several 
measurements for each frequency range were reported. The level of the electromagnetic 
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noise decreases by 20 dB for each decade of distance; thus, that level is less than the 
upper limit determined by technical requirements 6 when distances are greater than 100 m . 
Thus, the calculations show that the field strength depends on arc peak current, the 
current pulse form, and the distance from a discharge site. Because the arcing rate 
depends strongly on both bias voltage and plasma density 2 , an increase in operating 
voltage for solar cells will inevitably generate more intensive electromagnetic noise around 
the spacecraft, particularly in LEO conditions. In spite of the fact that the physical 
mechanism of arcing is not finally established yet, simple estimates can be made 



Fig. 6. Contour map of the electromagnetic field strength on log/-logx plane. 


/ = — r 1 = C • — 7 " , where Q is electrical charge, and V is the voltage drop between the 


regarding the expected peak current value. Really, an arc is an electrical disharge between 
two electrically charged bodies that can be considered as a capacitor with 
the capacitance C. 

The transient current is determined by the obvious relation 

f iQ_ r 

dt dt 
bodies. 

dV V 

The temporal rate of the voltage change can be estimated as — — » - J2!5 - . According to 

dt t 

SAMPIE experimental data the peak voltage (it is equal to the bias voltage) Kuu =400- 
600 V, and the average discharge time z=0. 1-0.5 fis . Now it is possible to estimate the 

C 

peak current / mmx = (0.8 - 6) • — A. The last result is in agreement with 

measured values (SAMPIE data, 1994). 
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Before beginning the computations of plasma waves it seems reasonable to show some 
results of the measurements of the arc peak current and duration of pulses that were 
obtained by SAMPIE and the LP experiment. The results of two different experiments 
(E 60-1 and E_64-l) are shown in Fig. 7 and Fig. 8 respectively. It is seen that the peak 
current was as high as I max =3 A, and the duration of pulse is equal to 0.2-0.4 ps for most 
cases. The electrical charge that was lost by the sample of solar cells depends on the 
duration of pulse, and it ranged usually between 0.2 and 0.4 pCl. 

It should be mentioned that arcs with very different parameters were observed. For 
example, during the experiment E48-1 five arcs were registered, and the peak current 
was no more than 0.6 A with approximately equal duration for each arc (0.7 ps). Only 
one arc was observed during the experiment E_44-l: the peak current was equal to 1.73 
A, and the duration was equal to 0.44 ps. Thus, the stochastic nature of the arcing onset 
demands consideration of the generated electromagnetic radiation within the framework 
of a statistical theory. The calculations above are useful as a first step to elaborate the 
theory of EMR generated by arcing on spacecraft surfaces. 



PULSE DURATION (MICROSEC) 


Fig.7a. Data from the Experiment #60-1. The bias voltage V— 400 V, and 181 arcs 
were registered. 
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Fig. 7b. Lost charge vs. pulse duration for the same experiment. A trend is 
clear (correlation coefficient corr(I, T)=0.259). 

In order to obtain more details regarding the arcing process and plasma wave 
characteristics ground experiments were carried out. 11 A plate of anodized aluminum 
(18x18 cm) was biased at -275 V and mounted inside a vacuum tank with neutral gas 
pressure p»20 (iTorr and electron number density n,~l(f cm* 3 . The electron temperature 
was approximately equal to 1 eV, and the ion temperature could be estimated as 300 K. 
Three Langmuir probes (A, B, and C) were placed in this tank at distances 55, 75, and 100 
cm from the plate respectively. Arc current and the potential of one of the Langmuir 
probes were measured simultaneously. Experimental data were represented in a digital 
form with a sampling interval of 5 ns that allowed us to analyze fluctuations with 
frequencies up to 100 MHz. 

As it was mentioned above, LP ground experiments were performed with an additional 
capacitor installed between a conducting plate and ground. The results of measurements 
are shown in Table 1. One can see from this table that there are no correlations between 
the electron number density and the value of electrical charge lost due to arcing. This 
means that we observed plasma waves in the background plasma, although there is some 
evidence of changes in electron number density in the vicinity of the Langmuir probes. 
We discuss these details below. The magnitude of the peak current increases with 
increasing of capacitance. This fact can be explained by the following estimates: peak 
current Ima^CV/At, and the discharge time At depends on plasma parameters and 
conducting plate material characteristics only. We should stress that the last statement is 
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valid for an arc current high enough to initiate the secondary processes in the arc column. 
For the initial stages of the discharge the rise time grows with increasing capacitance. 
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PULSE DURATION (MICROSEC) 


Fig. 8a. Data from the Experiment #64-1. The bias voltage V=-400 V, and 750 arcs 
were registered. 

6. High-frequency longitudinal waves 


The equation for the Fourier transformation of the electrical field strength can be written 
within the framework of a high frequency approximation (/»Fh , fw ) in the following 
form: 


Z / 2 P(f) b i -b J + S(f) • (£„. -b r bj)-i’ D(f) • • b k E,(f, k„k 2> k 3 ) - k 2 • £,(/, k , , k 2 , k 2 ) 

j =\ L *=i J 


+ • Z Kj • Ej (/ , K \ » ^2 » ^3 ) = - ] 1- 8 • 1 1 0 4 • ' i • f • ®( /) • a • . w, 

>=1 


( 17 ) 


where 


s if) = Ji /p2 ; = yi K 2 =c 2 ‘k 2 \ K = (k ] ,k 2 ,k 3 ). 
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and the electric field strength and frequency are expressed in the units Vm 3 /m*MHz and 
MHz respectively. 



Fig.8b Lost charge vs. pulse duration for the same experiment. A trend is 
clear (correlation coefficient corr(I, T)=0.42 7). 

The computations below are performed for the geometry shown in Fig. 9. 



Fig. 9 The geometry that was used for computations described below. 
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Table 1 









exper 

c. 

I max. 

At, 

finax. 

n./10 6 

Umax, 

TpulM 

AQ, 


(pF) 

(A) 

(ps) 

(MHz) 

(cm 3 ) 

(V) 

(ps) 

(PCI) 

LPA4 

0.15 

3.5 

9 

18 

4 

0.08 

7 

31.5 

LPA5 

0.15 

3 

10 

18 

4 

0.02 

5 

30 

LPB5 

0.15 

1.2 

15 

15 

3 

0.3 

10 

18 

LPA7 

0.15 

3 

10 



0.03 

10 

30 

LPB6 

0.15 

2.5 

12 

12 

1.5 

0.12 

7 

30 

LPB7 

0.15 

2.5 

15 

10 

1.2 

0.08 

10 

37.5 

LPA8 

0.15 

3 

10 

20 

5 

0.03 

10 

30 

LPB8 

0.15 

2.7 

10 

15 

3 

0.06 

10 

27 

LPA9 

0.15 

3 

12 

18 

4 

0.03 

7 

36 

LPB9 

0.15 

2 

13 

10 

1.2 

0.2 

8 

26 

LPB10 

0.15 

2.8 

10 

14 

2.4 

0.1 

7 

28 

DLPA2 

0.5 

12 

10 

17 

3.5 

0.2 

7 

120 

DLPA3 

0.5 

7 

17 

16 

3 

0.22 

8 

119 

DLPA8 

0.5 

13 

8 

20 

5 

0.05 

5 

104 

DLPB8 

0.5 

14 

8 

15 

3 

0.3 

7 

112 

DLPB2 

0.5 

8.2 

15 

12 

1.6 

0.75 

15 

123 

DLPB4 

0.5 

9 

13 

11 

1.4 

0.5 

14 

117 

DLPC8 

0.5 

9 

13 

12 

1.6 

0.15 

10 

117 

DLPC4 

0.5 

13 

8 

6 

0.7 

0.03 

10 

104 

DLPC2 

0.5 

8.2 

13 

14 

2.4 

0.15 

10 

107 

ELPA2 

0.22 

2.5 

12 

15 

3 

0.05 

7 

30 

ELPA3 

0.22 

2.8 

12 

15 

3 

0.06 

8 

33.6 

ELPA4 

0.22 

3.6 

10 

22 

6 

0.05 

8 

36 

ELPA5 

0.22 

4.1 

10 

17 

5 

0.18 

5 

41 

ELPA6 

0.22 

3.8 

9 

10 

1.2 

0.02 

10 

34.2 

ELPB1 

0.22 

3.3 

9.5 

9 

1 

0.3 

8 

31.4 

ELPB4 

0.22 

2.5 

14 

10 

1.2 

0.35 

10 

35 

ELPB5 

0.22 

4.3 

9.1 

9 

1 

0.75 

30 

39.1 

ELPC1 

0.22 

3.1 

13.6 

5 

0.2 

0.02 

5 

42.2 

ELPC2 

0.22 

2.6 

12 

8.8 

1 

0.45 

10 

31.2 

ELPC3 

0.22 

2.2 

16 

9 

1 

0.6 

13 

35.2 

ELPC4 

0.22 

2 

16 

13 

2 

0.075 

12 

32 

ELPC5 

0.22 

4.1 

10 

9 

1 

0.02 

10 

41 

ELPC6 

0.22 

2.3 

13 

12 

1.6 

0.1 

8 

30 

corr(Imax,Umax ) 

corr(C,AQ) 

corr(Umax,AQ) 

corr(n., AQ) 

corr(At,T pu i K1 ) 

0.1126199 


0.9791811 

0.2534364 

-0.001009 

0.096377 


B -probe 







I max At 

n. 

Umax 

TpulM 

AQ 

A(f mKt ) 




1.2 15 

3 

0.3 

10 

18 

72 




2.5 12 

1.5 

0.12 

7 

30 

31 




2.5 15 

1.2 

0.08 

10 

37.5 

36 




2.7 10 

3 

0.06 

10 

27 

54.5 




2 13 

1.2 

0.2 

8 

26 

62 




2.8 10 

2.4 

0.1 

7 

28 

55 




14 8 

3 

0.3 

7 

112 

26 




8.2 15 

1.6 

0.75 

15 

123 

130 




9 13 

1.4 

0.5 

14 

117 

55 




3.3 9.5 

1 

0.3 

8 

31.4 

20 




2.5 14 

1.2 

0.35 

10 

35 

25 




4.3 9.1 
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The correlational characteristics do not change significantly if we take into account data 
from one probe. This is seen in the row of numbers that are calculated for B probe only: 
corr(Imax,Umax) corr(At,T |)Ul *) corr(n.,A Q) corr(L.„.At£~.tt corrfUn^, A(f m<a )) 

0.414447 - 0.096776 0.0405188 0.0921994 0.5849428 



f 

Frequency (MHz) 

Fig.lOa The electric field strength derived from Eq.17. Hie resonance at the 

plasma frequency (f e =18 MHz ) is very strong for long longitudinal waves 
(kfI) travelling along the magnetic field. 

(iv=0.5 MHz , 9 =45° , and the attenuation rate y If =0.03) 


Now we have theoretical results that demonstrate the general features of plasma waves 
generated by arcing on the conductive plate immersed in a low density plasma, and we can 
compare experimental data with theoretical predictions. The results of one experiment are 
shown in Fig. 11. 

We may conclude that there is a quite good agreement between theoretical calculations 
(Fig. 10) and experimental data (Fig. lid) concerning the spectrum of plasma waves. 

To understand the origin of the second pulse in Fig. 1 lc we computed the spectrogram of 
fluctuations for this paricular experiment (Fig. 12). One can see that the frequencies of 
plasma oscillations during the first pulse lie in the range 10-12 MHz , and these 
frequencies increase up to 22 MHz during the second pulse (80-100 ps). If we suggest 
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Fig. 10b. The electric Held strength derived from Eq. 17 for longitudinal plasma 
waves traveling perpendicularly to the magnetic field. 

(F.=0.5MHz, <9=45°, k=5, and y/f t =0.03). 

The splitting of the maximum at the plasma frequency is seen 
very clearly. 


that Langmuir waves were registered in this experiment, the changes in the frequencies 
of oscillation should mean a considerable increase of the electron number density - from 
(1. 2-1.8)10* up to 6- 10* cm' 3 . This growth could be explained by the expanding arc 
plasma front that intersected the Langmuir probe. 18 Moreover, because the time interval 
between the two pulses is approximately equal to 25 ps (Fig. 11c), and the distance 
between LP and arc site is equal to L= 75 cm, we can calculate the speed of the plasma 


front , = — — « 30 km / s . This estimate for front velocity derived from our 
} Ar 


experimental data is in good agreement with the direct measurements 18 performed in a 
plasma with almost the same parameters as in our experiment. But the absence of 
correlations between the frequencies of oscillations and the magnitudes of lost electrical 
charges that were determined for the second pulses in a set of experiments does not allow 
us to be absolutely sure of the adequacy of this simple model. 
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Fig. 11a. Experiment LPB_9. The additional capacitance is equal to 0.15 pF, 
and the amplitude of current is equal to 2 A. 
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Fig. lib. The spectrum of the arc current pulse shown above. 
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Fig. 11c. The Langmuir probe potential fluctuations. Two separated pulses are 
seen quite clear. The first pulse has the highest amplitude, and its 
maximum is simultaneous to the maximum of the arc current pulse. 


The amplitude of the electron number density fluctuations can be estimated from the 
obvious relation: 

n = 0.02 -0.75 

«. 1 (18) 

To obtain the estimate (18) we used data from the Table 1 for Umax and measurements of 
the electron temperature which indicated the magnitude T e & 1 eV. It is much more 
difficult to obtain information regarding the wave numbers of observed waves. Really, 
the LP provides measurements that represent the fluctuations of electrostatic potential 
integrated over wavenumbers: 

U(f) = \u(f,k)d'k (19) 


21 


10 



Fig.lld. Spectrum of fluctuations shown above. Resonance at plasma frequency 
can be easily observed. 


Thus, we may determine the range of wave numbers only but not the spectral density of 
fluctuations U(f,k) itself The upper limit can be estimated from the formula for Landau 
damping of plasma waves. According to Ref.(13) the ratio of decrement to the plasma 
frequency can be written in the following form: 


- 7 - = G(x) = 2n • 




1 

2-x 1 


( 20 ) 


where x = k-A t . 

It seems convenient for further analysis to represent Eq. (20) in the graph (Fig. 13). If 
we suggest that Landau damping is the only mechanism for the attenuation of plasma 
waves, we can estimate the decrement as G(x )£ 0 .01 from the data shown in Fig. 1 lc. 
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This means that x<0.2 (Fig. 13), and the lower limit on the wave length = 10 • n • A, 



Fig. 13 . Attenuation rate vs. wave number for plasma waves. 

For example, during the experiment LPB_9 the electron number density was equal to 
n, = 10‘ cm 1 . The Debye length can be easily calculated : A, = 0.6 cm ; thus, we 
get the lower limit on the length of plasma waves ^ = 18 cm The phase speed of 
electrostatic waves V p = X- f exceeds 210 8 cm/s ; thus, we did not have any chance 

to register a time delay in the pulse arrival time even if we were performing measurements 
for two or three probes simultaneously. 

As it is seen from Table 1, the amplitude of fluctuations was very high for several 
experiments: t/ ma * > 0.3 V. We may expect that nonlinear conversion of plasma waves 
influences the time behaviour of pulses in such cases. One of the possible nonlinear 
processes is the generation of ion acoustic waves by the decay instability. 20 This process 
has the threshold - the critical magnitude of the electric field strength, shown in Fig. 14 
for different electron number densities and electron temperatures. 
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Fig. 14. Critical electric field strength vs. dimensionless wavenumber for two 
electron number densities #t e =10 n and 10 12 nf 3 and for two electron 
temperatures 7V=10 4 and 210 4 K. 
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Fig. 15 . Increment of plasma waves damping due to nonlinear conversion for three 
magnitudes of potential ij= 0.1, 0.5 and 1. 


One can see in the graph above that the increment of nonlinear conversion becomes as 
large as the increment of Landau damping for waves long enough (x<0.2) when the 
amplitude of fluctuations exceeds ^ cr =0.5. However, ion acoustic waves could not be 
registered even during three experiments with the largest amplitudes of fluctuation 
(DLPB2, DLPB_4, and ELPB_5) because of the low frequency of such waves (a few 
kHz) . We might suggest that some nonlinear processes took place in the plasma during 
the experiments with the highest level of fluctuations. Such nonlinear processes should 
cause a widening of the resonance curve at the plasma frequency and the generation of 
additional broadband noise, and these expectations are confirmed by measurements (see 
Fig. 16). On the contrary, experiments with low level of fluctuations demonstrate very 
narrow resonance at the plasma frequency (one example is shown in Fig. 17). 

Now we can estimate the amplitude of electric field strength for longitudinal plasma 
waves: 


£=*•£/_ 


2 -x-U 


( 22 ) 


During the experiment ELPB5 , Umax-' 0.75 V and /l m/n =0.2 m; thus, we obtain an upper 
limit on the electric field strength iw=24 V/m. This is a very rough estimate because 
each experiment demonstrates a complicated spectral distribution of fluctuations. For the 
experiment DLPA 3, we estimate the spectral density of electric field as E(f)~ 1 
V/m-MHz in the vicinity of the resonance (Fig. 17). 
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Fig. 16. Experiment ELPB_5. Spectrum of fluctuations differs significantly from 
one in Fig. lid perhaps because of nonlinear processes developed in 
the plasma with a high level of fluctuations. 


7. Conclusion 


One of the most important problems facing large spacecraft with powerful high voltage 
solar arrays is the generation of electromagnetic radiation by arcing on its surfaces. As it 
follows from our calculations and available experimental data, broadband electromagnetic 
noise with electric field amplitudes as high as a few volts per meter can be generated by 
arcing. To prevent a negative influence of this radiation on spacecraft operation additional 
work should be done: i)to establish technical requirements on equipment and devices that 
could be influenced; ii)to determine the spectrum and level of radiation that could be 
expected in the vicinity of the equipment; iii)to elaborate methods for mitigating adverse 
consequences. 
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Fig. 17. Experiment DLPA_3. Amplitude of fluctuations is relatively low ( 7 = 0 . 2 ), 
and spectrum demonstrates a narrow resonance at the plasma frequency. 
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